Unification of perturbation theory, RMT and semiclassical considerations in the study 

of parametrically-dependent eigenstates 
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We consider a classically chaotic system that is described by an Hamiltonian Tt(Q,P;x) where 
a; is a constant parameter. Our main interest is in the case of a gas-particle inside a cavity, where 
x controls a deformation of the boundary or the position of a 'piston'. The quantum-eigenstates of 
the system are \n(x)). We describe how the parametric kernel P(n\m) = |(n(a;)|m(xo))| 2 evolves 
as a function of Sx — x—xq. We explore both the perturbative and the non-perturbative regimes, 
and discuss the capabilities and the limitations of semiclassical as well as of random-waves and 
random-matrix-theory (RMT) considerations. 



Consider a system that is described by an Hamiltonian 
Ti(Q, P; x) where (Q, P) are canonical variables and x is 
a constant parameter. Our main interest is in the case 
where the parameter x represents the position of a small 
rigid body ('piston') which is located inside a cavity, and 
the (Q, P) variables describe the motion of a 'gas parti- 
cle'. It is assumed that the system is classically chaotic. 
The eigenstates of the quantized Hamiltonian are \n(x)) 
and the corresponding eigen-energies are E n (x). The 
eigen-energies are assumed to be ordered, and the mean 
level spacing will be denoted by A. We are interested in 
the parametric kernel 

|(ra(a;)|m(xo))| 2 = trace(p„p m ) 



P(n\m) 



(1) 



In the equation above p m (Q,P) and p n (Q,P) are the 
Wigner functions that correspond to the eigenstates 
\m(xo)) and \n(x)) respectively. The trace stands for 
dQdP/ (2irh) d integration. The difference x — xq will be 
denoted by Sx. We assume a dense spectrum. The ker- 
nel P(n\m), regarded as a function of n— to, describes 
an energy distribution. As Sx becomes larger, the width 
as well as the whole profile of this distribution 'evolves'. 
Our aim is to study this parametric-evolution (PE). 

The understanding of PE is essential for the analy- 
sis of experimental circumstances where the 'sudden ap- 
proximation' applies fm. It also constitutes a prelimi- 
nary stage in the studies of quantum dissipation Jt| . The 
function P(n\m) has received different names such as 
'strength function' j| and 'local density of states' |Q. 
Some generic features of PE can be deduced by refer- 
ring to time-independent first-order perturbation theory 
(FOPT), and to random-matrix-theory (RMT) consider- 
ations |^||] . Other features can be deduced using classical 
approximation ^J,^|, or its more controlled version that 
we are going to call phase-space semiclassical approxi- 
mation 0. Still another strategy is to use time-domain 
semiclassical considerations Q|. In case of cavities one 
can be tempted to use 'random-wave' considerations as 
well. Depending on the chosen strategy, different results 
can be obtained. The 'cavity' system is a prototype ex- 
ample for demonstrating the 'clash' between the various 
approaches to the problem. 

We are considering the cavity example where we have 
a 'gas' particle whose kinetic energy is E — ^mv 2 , where 



to is its mass, and v is its velocity. The 'gas' particle is 
moving inside a cavity whose volume is V and whose di- 
mensionality is d. The ballistic mean free path The 
area of the displaced wall-element ('piston' for brevity) 
is A, while its effective area is A^t, see M for geometrical 
definition. The mean free path £ cal ps \l j A between col- 
lisions with the piston may be much smaller compared 
with t bl . The penetration distance upon a collision is 
I = E/ /, where / is the force that is exerted by the wall. 
Upon quantization we have an additional length scale, 
which is the De-Broglie wavelength A B = 2ttTi / (mv) . We 
shall distinguish between the hard walls case where we 
assume £ < A B <C £ bl , and soft walls for which A B <C £■ 
Note that taking Ti — * implies soft walls. 

For convenience of the reader we start by listing the 
various expressions that can be derived for P(7i|to), along 
with an overview of our PE picture. Then we proceed 
with a detailed presentation. We are going to argue that 
four parametric scales Sx^ <§; <5x NU <C 5x VIt Sx sc are 
important in the the study of PE. 

Standard FOPT assumes that P(n\m) has a simple 
perturbative structure that contains mainly one state: 



P(?i|to) ps 5 nm + Tail(n — to) 



(2) 



We define Sx^" 1 to be the parametric change that is re- 
quired in order to mix neighboring levels. For Sx > Sx^ 
an improved version of FOPT implies that P(n|r 
core-tail structure |tJ : 



has a 



P(n|m) w Core(n — m) + Tail(rt — to) 



(3) 



The core consists of those levels that are mixed non- 
perturbatively, and the tail evolves as if standard FOPT 
is still applicable. In particular we argue that the tail 
grows like Sx 2 , and not like Sx. We also explain how 
the core-width depends on Sx. It should be noted that 
Wigner's Lorentzian can be regarded as a special 
case of core-tail structure. 

Another strategy is to use semiclassical considerations. 
The simplest idea is to look on the definition (Q) and to 
argue that p n (Q,P) and p m (Q,P) can be approximated 
by microcanonical distributions. This is equivalent to 
the classical approximation that has been tested in ||. 



1 



If we try to apply this approximation to the cavity ex- 
ample we should be aware of a certain complication that 
is illustrated in Fig.l. One obtains 



P(n\m ) 



S(n — m) + S 



E n — En 
SE cl 



(4) 



The detailed explanation of this expression is postponed 
to a later paragraph. A more careful semiclassical proce- 
dure is to take the width of Wigner function into account. 
Namely, we can approximate p n (Q,P) and p m (Q,P) by 
smeared microcanonical distributions. It can be used in 
order to get an idea concerning the quantum mechanical 
'interpretation' of the Dirac's delta function component 
in (El). The result is 



S(n — m) 



5E S 



1 



7T SE 2 C + {E n ~E m ) 2 



(5) 



with <5-E sc = h/r h i, where r bl = l bl /v. However, we are 
going to argue that the latter procedure, which is equiv- 
alent to the assumption of having uncorrelated random 
waves, is an over-simplification. It is better to use the 
time-domain semiclassical approach which is based on 
the realization that P{n\m) is related to the so-called 
survival amplitude via a Fourier transform JjJ , leading to 
the identification SE SC = ?i/r col , where r col = £ col /u. 
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FIG. 1. Phase space illustration of the energy surface (represented 
by bold solid line) that support the Wigner function of a given 
eigenstate \m(xo)), and the energy surfaces (light solid lines) that 
support the Wigner functions of some of the eigenstates \n(x)). 
The left illustration refers to an hypothetical generic case, while 
the right illustration refers to the cavity example. The associated 
P(n\m) is plotted below each of the phase space illustrations: The 
classical behavior is indicate by the black lines, and the quantum- 
mechanical behavior is represented by the grey filling. It should 
be realized that detailed quantal-classical correspondence assumed, 
which is guaranteed only if 8x > Sxsc- In the quantum-mechanical 
case classical sharp-cutoffs are being smeared ('tunneling' correc- 
tion). In the cavity example the classical delta-singularity is being 
smeared as well. In the latter case a naive phase-space picture 
cannot be used in order to determine the width of the smearing. 



The important point to realize is that (0) with (^|) is 
fundamentally different from either (ph or . The main 
purpose of this Letter is to give a clear idea of the route 
from the regime where perturbation theory applies to the 
non-perturbative regime regime where semiclassical con- 
sideration become useful. We are going to explain that 
the width of the core in (|^) defines a 'window' through 
which we can view the 'landscape' of the semiclassical 
analysis. As Sx becomes larger, this 'window' becomes 
wider, and eventually some of semiclassical structure is 
exposed. This is marked by the non-universal paramet- 
ric scale 5x NV . For Sx much larger than S j the non- 
universal structure (ph of the core is exposed. Still, the 
perturbative tail of (y) may survive for relatively large 
values of Sx. One wonders whether this tail survives 
for arbitrarily large Sx. While the answer for the latter 
question may be positive for hard walls, it is definitely 
negative for soft walls, as well as for any other generic 
system. Assuming soft walls, one should realize that the 
perturbative tail of (||) occupies a finite bandwidth. It is 
well known || that having finite bandwidth is a generic 
feature of all quantized systems, provided Ti is reason- 
ably small. Therefore one should introduce an additional 
parametric scale Sx plt . For Sx 3> Sx pvt the core spills 
over the bandwidth of the perturbative tail, and P(n\m) 
becomes purely non-perturbative. The non-perturbative 
P(n\m) does not necessarily correspond to the classical 
approximation (EJ). We are going to introduce one more 
additional scale Sx sc - For Sx 3> Sx sc detailed quantal- 
classical correspondence is guaranteed, and (|J) with (||) 
becomes applicable. 

Expression (||) is a straightforward result of standard 
time-independent FOPT where 



Tail(n 



dU 

dx 



Sx 2 



(E n —E m ) 2 



(6) 



An estimate for the matrix elements {dTL/ dx) nm follows 
from simple considerations (?[. Upon substitution into 
(p) it leads to: 



P(n\m) 



( Sx 



1 



ySxT 1 J (n — m) 2+1 
for b(x) <C \n — m\ <C b 



with /3 = 2 and b(x) — 0. We have defined 



'r((d+3)/2) i Ad+1 



4 7r (d-i)/2 a 



(7) 



(8) 



We shall refer to the dependence of \(d7i/ dx) nm \ 2 on 
n — m as the band-profile. It is well known that the 
band-profile is related (via a Fourier transform) to a clas- 
sical correlation function. If successive collisions with the 
'piston' are uncorrelated then we have 7=0. But in other 
typical circumstances || we may have 0<7. The matrix 
(dH/dx) nm is not a banded matrix unless we assume soft 
(rather than hard) walls. In the latter case the bandwidth 
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Ab = (h/T cl ) is related to the collision time r cj = IJv with 
the walls. Having hard walls (I < A B ), implies that A;, 
becomes (formally) larger than E. The notion of band- 
width is meaningful only for soft walls (I 3> A B ). In di- 
mensionless units the bandwidth it is commonly denoted 
by b = Aft/A. 

The standard result (|) with (g) of FOPT is valid as 
long as Sx <C fag m . Once Sx becomes of the order of Sx*" 1 , 
we expect few levels to be mixed non-perturbatively. 
Consequently (for Sx > Sx*" 1 ) the standard version of 
FOPT breaks down. As 5x becomes larger, more and 
more levels are being mixed non-perturbatively, and it is 
natural to distinguish between core and tail regions. The 
core-width b(x) is conveniently defined as the participa- 
tion ratio (PRR), namely ^MEnC^I™)) 2 ) -1 . Thc 
tail consists of all the levels that become 'occupied' due 
to first-order transitions from the core. It extends within 
the range b(x)<\n— m\<b. Most of the spreading proba- 
bility is contained within the core region, which implies 
a natural extension of FOPT: The first step is to make 
a transformation to a new basis where transitions within 
the core are eliminated; The second step is to use FOPT 
(in the new basis) in order to analyze the core-to-tail 
transitions. Details of this procedure are discussed in [7\, 
and the consequences have been tested numerically 
The most important (and non-trivial) consequence of this 
procedure is the observation that mixing on small scales 
does not affect the transitions on large-scales. Therefore 
we have in the tail region P(n\m) oc Sx 2 rather than 
P(n\m) oc Sx. The above considerations can be summa- 
rized by stating that ^ holds with /3 = 2 well beyond 
the breakdown of the standard FOPT. 

We turn now to discuss the non-perturbative structure 
of the core. The identification of b(x) with the inverse- 
participation-ratio is a practical procedure as long as we 
assume a simple energy spreading profile where the core 
is characterized by a single width-scale. As long as this 
assumption (of having structure-less core) is true we can 
make one step further and argue that 



b(x) 



2n 2 



Sx 
SxT' 



2/(1+7) 



assuming \j\ < 1 (9) 



The argument goes as follows: Assuming that there 
is only one relevant energy scale (b(x)) it is im- 
plied by (R) that P(n\m) has the normalization 
(Sx I 'Sx^" 1 ) 2 1 '(b(x)) 1+ ' f . This should be of order unity. 
Hence (||) follows. The tail should go down fast enough 
(7 > —1) else our 'improved' perturbation theory does 
not hold. Namely, for 7 < —1 the core- width becomes 
cutoff dependent (via its definition as an PRR) , and con- 
sequently it is not legitimate to neglect the 'back reac- 
tion' for core-to-tail transitions. The tail should go down 
slow enough (7 < 1) in order to guarantee that the core 
width is tail-determined. Else, if 7 > 1 then the core 
width is expected to be determined by transitions be- 
tween near-neighbor levels leading to a simple linear be- 
havior b(x) = (Sx/Sx*" 1 ). 



Non-perturbative features of P(n\m) are associated 
with the structure of the core. In order to further analyze 
the non-perturbative features of P(n\m) we are going to 
apply semiclassical considerations. An eigenstate \n(x)) 
can be represented by a Wigner function p n (Q,P). In the 
classical limit p n (Q,P) is supported by the energy sur- 
face ri(Q,P;x) = E n . However, unlike microcanonical 
distribution, it is further characterized by a non-trivial 
transverse structure. One should distinguish between the 
'bulk' flat-portions of the energy-surface (where Q de- 
scribes free motion), and the relatively narrow curved- 
portions (where Q is within the wall field-of- force) . In 
the curved-portion of the energy surface (near the turn- 
ing points), Wigner function has a transverse Airy struc- 
ture whose 'thickness' is characterized by the energy scale 
A sc = ((H/t^E) 1 ' 3 . This latter expression is valid for 
soft walls (A B C In the hard wall case (I < A B ) it 
goes to A sc ~ E. Unlike the curved-portions, the 'bulk' 
flat-portions of the energy surface are characterized by 
A sc = (h/i~ bl ). Now we consider two sets of eigenstates, 
\n(x)) and \m(xo)), which are represented by two sets of 
Wigner functions p n (Q, P) and p m (Q, P). The probabil- 
ity kernel (|l|) can be written as P(n\m) — trace(p„p m ). 
If Pn{Q,P) and p m (Q,P) are approximated by micro- 
canonical distributions, then P(n\m) is just the projec- 
tion of the energy surface that correspond to m, on the 
"new" energy surface that correspond to n. This leads to 
the classical approximation Eq.(^). In the classical limit 
n and m become continuous variables, and Dirac's delta 
just reflects the observation that most of the energy sur- 
face (the 'bulk' component) is not affected by changing 
the position of the classically-small 'piston'. The second 
term in (||) has the normalization (t cJ /t coJ ), and corre- 
sponds to the tiny component which is affected by the 
displacement of the 'piston'. For S x < I it extends over 
an energy range SE cl = fSx, where / is the force which 
is exerted on the particle by the wall. When Sx becomes 
larger than t the energy spread becomes of order E. 

In the quantum-mechanical case we should wonder 
whether (Q) can be used as an approximation, and what 
is the proper 'interpretation' of Dirac's delta function. 
It is relatively easy to specify sufficient condition for 
the validity of the classical approximation. Namely, the 
transverse structure of Wigner function can be ignored if 
A sc <C \E n — E m \ -C 5E cl . For hard walls A sc ~ E and 
therefore the classical approximation becomes inapplica- 
ble. For soft walls the necessary condition A sc <c 8E cl is 
satisfied provided Sx is large enough. Namely Sx ^> Sx sc , 

with sx sc = 

We want to go beyond the classical approximation, 
and to understand how the classical Dirac's delta func- 
tion in (fj) manifests itself in the quantum mechanical 
case. Thus we are interested in the singular overlap 
of the 'bulk' components (see Fig.l), and the relevant 
A sc for the current discussion is h/r bl . The most naive 
guess is that the contribution due to the overlap of 'bulk' 
components becomes non-zero once \E n — E m \ < A sc . 
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Equivalently, one may invoke a 'random-wave' assump- 
tion: One may have the idea that \n(x)) and |m(a;o)) 
can be treated as uncorrelated random-superpositions of 
plane-waves. Adopting the random-wave assumption, it 
is technically lengthy but still straightforward to derive 
(@) with SE SC = U/t m . 

The naive phase-space argument that supports the 
'random wave' result (|^) is definitely wrong. One should 
realize that \E n — E m \ < A sc is a necessary rather than 
a sufficient condition for having a non-vanishing 'bulk' 
contribution. This latter observation becomes evident if 
one considers the trivial case Sx = for which we should 
get P(n\m) = for any n ^ m. Thus h/r bl should be 
regarded as an upper limit for SE SC . We are going to 
argue that the correct result (for large enough Sx) is in- 
deed (§), but T b i should be replaced by the possibly much 
larger length-scale r col . 

In order to go beyond the random-wave assumption 
we use the time-domain semiclassical approach which is 
based on the realization that P(n\m) is related to the 
so-called survival amplitude via a Fourier transform 0]: 

P(n\m)2ir5(u>-^) = FT (m\ exp(-^)|m) (10) 

n 

Note that \m) is an eigenstate of H(Q, P; xq) while 
H = H(Q,P;x). The knowledge of the short time dy- 
namics, via classical considerations, can be used to obtain 
the 'envelope' of P(n\m). Adopting Wigner's picture, 
the evolving |m) in the right hand side of ([n]) is rep- 
resented by an evolving (quasi) distribution p m (Q, P;t). 
Let us assume that the 'piston' is small, such that the 
collision rate with it (1/t co1 ) is much smaller than 1/t m . 
Due to the chaotic nature of the motion successive colli- 
sions with the piston are uncorrelated. It follows that the 
portion of p m {Q, P;t) which is not affected by collisions 
with the 'piston' decays exponentially as exp(— t/r col ). 
It is reasonable to assume that any scattered portion 
of p m (Q, P;t) lose phase-correlation with the unscat- 
tered portion. Therefore the right hand side of ((!(]) is 
the Fourier transform of an exponential. Consequently 
P(n\m) should have the Lorentzian shape (|J), but the 
correct energy- width is SE SC = n/r col rather than U/r bl . 

For an extremely small parametric change such that 
Sx <C Sx^," 1 we have the simple perturbative structure 
(||). Then, for larger values of Sx the energy distribution 
develops a core. As long as this core is structure-less it is 
characterized by the single width-scale b(x) of @. Now 
we would like to define a new parametric scale <5a; NU . By 
definition, for Sx 3> Sx NU non-universal features mani- 
fest themselves, and the core is characterized by more 
than one width-scale. For our 'cavity' example this hap- 
pens once the semiclassical Lorentzian structure (||) is 
exposed. This happens when b(x) of (^|) becomes larger 
than SE SC /A, leading to 

1/2 



1 

4tt 



(d+1) 



A 

^4 ff 



(11) 



regime Sx 3> Sx NV , where the semiclassical Lorentzian 
structure (||) is exposed. It is also important to realize 
that in the non-universal regime we do not have a theory 
for the b(x) of (^). The derivation of (|^) is based on the 
assumption of having a structure-less core, and therefore 
pertains only to the universal regime. 

It is well known j(| that for any quantized system 
{dTL/ dx) nm is characterized, for sufficiently small h~, by a 
finite bandwidth Afc. Consequently it is possible to define 
a non-perturbative regime Sx Sx prt , where the condi- 
tion b(x) <g; 6 is violated. In the non-perturbative regime 
expression (Q) becomes inapplicable because the core 
spills over the (perturbative) tail region. Thus P(n\m) 
becomes purely non-perturbative. Hard walls are non- 
generic as far as the above semiclassical considerations 
are concerned. In the proper classical limit all the clas- 
sical quantities should be held fixed (and finite), while 
making h smaller and smaller. Therefore the proper clas- 
sical limit implies soft walls (A B <C I), leading to finite 
bandwidth Aj, = Ti/t c1 . From (||) it follows that the con- 
dition b(x) <C b is definitely not violated for Sx < Sx NV . 
Hence we conclude that Sx plt ^> Sx NV , but we cannot 
give an explicit expression since @ becomes non- valid in 
the non-universal regime. 

In the parametric regime <5x NU <C Sx <C Sx plt we have 
on the one hand 8E cl 3> A;,, and on the other hand 
b{x) <C b by definition. Therefore we cannot get in 
this regime a contribution that corresponds to the sec- 
ond term in (Q) . A necessary condition for the manifes- 
tation of this second term is Sx 3> Sx VIt . However, it 
should be realized that having Sx 3> Sx prt is not a suffi- 
cient condition for having detailed correspondence with 
the classical approximation. For our 'cavity' example 
detailed correspondence means that the whole classical 
structure of @) is exposed. As discussed previously, the 
sufficient condition for having such detailed correspon- 
dence is Sx ^s> Sx sc . This latter condition is always sat- 
isfied in the limit h — > 0. 

We thank Michael Haggerty for his useful comments, 
and ITAMP for their support. 
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Let us re-emphasize that the semiclassical argument that 
is based on ( |l0| ) applies to the non-universal parametric 
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